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LISTING OF CLAIMS 

This listing of claims will replace all prior versions, and listings, of claims in the 
instant application. 

1 . (Currently amended) A method for controlling hourglass deformations of an 
under-integrated solid element in a finite element analysis for designing and 
analyzing a structural product, the method comprising: 

establishing a local initial element coordinate system of the under-integrated 
solid element fe ^from an initial undeformed geometry of the solid element; 

establishing a local current element coordinate system of the solid element fef 
from a current deformed geometry of the solid element; 

calculating a set of initial nodal coordinates of the solid element in the local 
initial element coordinate system; 

calculating a set of current nodal coordinates of the solid element in the local 
current element coordinate system; 

evaluating a set of hourglass shape vectors of the solid element from the 
initial nodal coordinates; 

calculating a set of hourglass deformation magnitudes of the solid element 
us i ng tho sot of hourglass shap e v e ctors, and d i fference betw ee n the i nit i a l noda l 
coordinates and th e current noda l coord i nates of corn e r nodes of the solid 
olom e n t as set forth in Equation as follows: 

N 

Equation: g ia = Y,y*J i*u ' x u ) 

7=1 

where: 

N = total number of the corner nodes in the solid element. 
J = one of the corner nodes. 
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g ia - one of the hourglass deformation magnitudes (spatial direction i for 
hourglass mode a). 

Y a j ~ one of the hourglass shape vectors (node J for hourglass mode a), 

X if = one of the initial nodal coordinates (node J in spatial direction i of the 
local initial element coordinate system). 

x u = one of the current nodal coordinates (node J in spatial direction i of the 

local current element coordinate system) ; 

evaluating a set of generalized hourglass forces from the hourglass 
deformation magnitudes, the local initial nodal coordinates, and material constants of 
the solid element; and 

calculating a set of counter nodal forces for controlling the hourglass 
deformations in the local current element coordinate system from the generalized 
hourglass forces and the hourglass shape vectors, wherein the set of counter nodal 
forces is applied in directions opposing to us e d to res i st the hourglass deformations 
such that the hourglass deformations are controlled in the finite element analysis of 
for designing and analyzing a t he-structural product. 

2. (Previously presented) The method as recited in claim 1, further comprising: 

transforming the set of counter nodal forces from the local current element 
coordinate system to global coordinate system before adding to global force array. 

3. (Original) The method as recited in claim 1, further comprising: 

calculating all terms of an element stabilization matrix for the solid element 
from the hourglass shape vectors, the initial nodal coordinates, and material 
constants of the solid element. 
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4. {Previously presented) The method as recited in claim 3, further comprising: 

transforming the stabilization matrix from the local initial element coordinate 
system to global coordinate system before adding the terms of the stabilization 
matrix into global stiffness matrix. 

5. (Original) The method as recited in claim 1 , wherein the solid element is chosen 
from the group consisting of three-dimensional 8-node hexahedral element, 6-node 
three-dimensional pentahedral element, two-dimensional 4-node plane strain 
element and two-dimensional 4-node axisymmetric continuum element. 

6. (Currently amended) A software product embodied in a tangible computer 
readable storage medium and executing in a computing device for controlling 
hourglass deformations of an under-integrated solid element in a finite element 
analysis for designing and analyzing a structural product , the software product 
comprising: 

program code for establishing a local initial element coordinate system of the 
under-integrated solid element fer -from an initial undeformed geometry of the solid 
element; 

program code for establishing a local current element coordinate system of 
the solid element fo^from a current deformed geometry of the solid element; 

program code for calculating a set of initial nodal coordinates of the solid 
element in the local initial element coordinate system; 

program code for calculating a set of current nodal coordinates of the solid 
element in the local current element coordinate system; 

program code for evaluating a set of hourglass shape vectors of the solid 
element from the initial nodal coordinates; 

program code for calculating a set of hourglass deformation magnitudes of 
the solid element as set forth in Eguation as follows: 
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N _ 

Equation: g ia = ]T y tXj (x u -X u ) 
where: 

N = total number of the corner nodes in the solid element, 
J = one of the corner nodes, 

g Ux - one of the hourglass deformation magnitudes (spatial direction i for 
hourglass mode a), 

v r , = one of the hourglass shape vectors (node J for hourglass mode a). 

X u - one of the initial nodal coordinates (node J in spatial direction i of the 
local initial element coordinate system). 

i f j = one of the current nodal coordinates (node J in spatial direction i of the 

local current element coordinate svstem) uc i ng tho sot of hourg l ass shapo vectors, 
and d i ff e renc e b e tw ee n th e in i t i a l nodal coord i nates and tho curr e nt noda l 
coord i nat e s of corn e r nodes of th e solid e l ement ; 

program code for evaluating a set of generalized hourglass forces from the 
hourglass deformation magnitudes, the initial nodal coordinates, and material 
constants of the solid element; and 

program code for calculating a set of counter nodal forces for controlling the 
hourglass deformations in the local current element coordinate system from the 
generalized hourglass forces and the hourglass shape vectors, wherein the set of 
counter nodal forces is applied in directions opposing to uood to rooiot the hourglass 
deformations such that the hourglass deformations are controlled in the finite 
element analysis e £for designing and analyzing a t he-structural product, 

7. (Currently amended) The software product as recited in claim 6, further 
comprising: 
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program code for transforming the set of counter nodal forces from the local 
current element coordinate system to global coordinate system before adding to 
global force array. 

8. (Original) The software product as recited in claim 6, further comprising: 

program code for calculating all terms of an element stabilization matrix for 
the solid element from the hourglass shape vectors, the initial nodal coordinates, 
and material constants of the solid element. 

9. (Previously presented) The software product as recited in claim 8, further 
comprising: 

program code for transforming the stabilization matrix from the local initial 
element coordinate system to global coordinate system before adding the terms of 
the stabilization matrix into global stiffness matrix. 

10. (Original) The software product as recited in claim 6, wherein the solid element is 
chosen from the group consisting of three-dimensional 8-node hexahedral element, 
6-node three-dimensional pentahedral element, two-dimensional 4-node plane strain 
element and two-dimensional 4-node axisymmetric continuum element. 

1 1 . (Currently amended) A system for controlling hourglass deformations of an 
under-integrated solid element in a finite element analysis for designing and 
analyzing a structural product, the system comprising: 

an I/O interface; 

a data communications interface; 

a memory for storing computer readable code for an application module; 
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at least one processor coupled to the memory, the I/O device and the data 
communications interface, said at least one processor executing the computer 
readable code in the memory to cause the application module to perform operations 
of: 

establishing a local initial element coordinate system of the under-integrated 
solid element fe M rom an initial undeformed geometry of the solid element; 

establishing a local current element coordinate system of the solid element fef 
from a current deformed geometry of the solid element; 

calculating a set of initial nodal coordinates of the solid element in the local 
initial element coordinate system; 

calculating a set of current nodal coordinates of the solid element in the local 
current element coordinate system; 

evaluating a set of hourglass shape vectors of the solid element from the 
initial nodal coordinates; 

calculating a set of hourglass deformation magnitudes of the solid element as 
set forth in Equation as follows: 

n _ 
Equation: g ia = ]T faj (x Lf - X u ) 
j=\ 

where: 

N = total number of the corner nodes in the solid element. 
J - one of the corner nodes, 

2 ia - one of the hourglass deformation magnitudes (spatial direction i for 
hourglass mode a). 

Y uJ - one of the hourglass shape vectors (node J for hourglass mode a), 

X u - one of the initial nodal coordinates (node J in spatial direction i of the 
local initial element coordinate system). 
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jr., = one of the current nodal coordinates (node J in spatial direction i of the 

local current element coordinate svstero) using tho sot of hourglass shape vectors, 
and difforonco botwoon the i n i t i al nodal coordinates and tho current nodal 
coord i nat e s of corn e r nod e s of the solid e le ment ; 

evaluating a set of generalized hourglass forces from the hourglass 
deformation magnitudes, the initial nodal coordinates, and material constants of the 
solid element; and 

calculating a set of counter nodal forces for controlling the hourglass 
deformations in the local current element coordinate system from the generalized 
hourglass forces and the hourglass shape vectors, wherein the set of counter nodal 
forces is applied in directions opposing to usod to resist the hourglass deformations 
such that the hourglass deformations are controlled in the finite element analysis «f 
for designing and analyzing a t ke-structural product 

12. (Previously presented) The system of claim 1 1 , further comprising operations of: 

transforming the set of counter nodal forces from the local current element 
coordinate system to global coordinate system before adding to global force array. 

13. (Previously presented) The system of claim 11, further comprising operations of: 

calculating all terms of an element stabilization matrix for the solid element 
from the hourglass shape vectors, the initial nodal coordinates, and material 
constants of the solid element. 

14. (Previously presented) The system of claim 13, further comprising operations of: 

transforming the stabilization matrix from the local initial element coordinate 
system to global coordinate system before adding the terms of the stabilization 
matrix into global stiffness matrix. 
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1 5. (Previously presented) The system of claim 1 1 , wherein the solid element is 
chosen from the group consisting of three-dimensional 8-node hexahedral element, 
6-node three-dimensional pentahedral element, two-dimensional 4-node plane strain 
element and two-dimensional 4-node axisymmetric continuum element- 
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